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Abstract 

The charge transport behaviors of parallel double single electron transistors (SETs) 
are investigated by the Anderson model with two impurity levels. The nonequilibrium 
Keldysh Green's technique is used to calculate the current-voltage characteristics of 
system. For SETs implemented by quantum dots (QDs) embedded into a thin SiC"2 
layer, the interdot Coulomb repulsion is more important than the interdot electron 
hopping as a result of high potential barrier height between QDs and Si02 barrier. 
We found that the interdot Coulomb repulsion not only leads to new resonant levels 
but also creates negative differential conductances. 
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I. INTRODUCTION 



The transport properties of a single electron transistor (SET) have been experimentally and theo- 
retically studied by many groups [1-5]. Moreover, the interesting phenomena of SETs such as Kondon 
effect and Coulomb blockade have been demonstrated experimentally [1-3]. SETs comparing with 
conventional transistors offer the advantages of good scalability, ultralow power operation for future 
memory and quantum information applications. Recently, some Si or Ge-dot SETs fabricated by 
CMOS-compatible process techniques have been proposed and large Coulomb-blockade oscillations at 
room temperature have also been reported [6-9]. 

For practical integrated-circuit applications, it is crucial to clarify that if the current-voltage 
characteristics of a single SET will be affected by the nearby SETs or quantum dots (QDs). The 
main purpose of this work is to theoretically study the charge transports of SETs sitting in parallel 
as depicted in Fig. 1, where A-SET (with dot A) is regarded as a detector SET and B-SET (with 
dot B) is as a target SET. Each SET consists of a single dot and two electrodes (source and drain). 
The mechanism of interactions between two SETs mainly arises from the interdot interactions such as 
interdot Coulomb repulsion and interdot electron hopping. The latter is negligible as a result of high 
Si0 2 barrier. Therefore, we take into account only the interdot Coulomb repulsion in the calculation 
of current-voltage characteristics. Kim and Hershfield studied the interdot spin-spin interaction effect 
on the tunnelling current for the parallel double QDs in the Kondo regime [10]. On the contrary, this 
interdot spin-spin interaction is negligible for the Coulomb blockade regime and therefore, it is not 
taken into account in this study. 

II. ENERGY LEVELS AND PARTICLE INTERACTIONS 

To construct the Anderson model simulating the system shown in Fig. 1, we calculate the energy 
levels, intradot and interdot Coulomb interactions, and interdot electron hopping strengths using the 
effective mass model. These physical parameters are important in the calculation of current-voltage 
characteristics. Two cylindrical Ge QDs are embedded into a Si0 2 layer with a finite width W. The 
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layer is then placed in contact with metals to form an n-i-n structure. The metallic electrodes are 
considered since they can avoid carriers frozen occurred in the electrodes at low temperature [11]. To 
calculate the electronic structures of QDs, we consider the cylindrical Ge QDs with a radius R and 
h = 2R . The conduction band of Ge has four equivalent valleys aligned along [1,1,1], [1,-1,-1], [-1,1,- 
1], and [-1,-1,1] [12]. In this calculation, the multi-valleys and image charge effects are ignored for 
simplicity. Within the effective mass model, the Hamiltonian is 

r - h 2 



2m* (p, z) 



V + V QDti (p,zM(T) (1) 



where m*(p, z) denotes the position-dependent electron effective mass. Hence, m* s = m e for Si0 2 
and ttiqI = l/3(l/m^ + 2/m t ) for Ge, where me = 1.59m e and m t = 0.0823m e are the electron 
longitudinal and transverse effective mass of bulk material. Vq_d,i(p, z) is approximated by a constant 
potential Vq in the QD region. Its value is determined by the conduction band offset between Ge and 
Si0 2 - Besides the energy levels, the Coulomb interactions should not be neglected for the electron 
transport through small QDs. The Coulomb interactions are calculated according to 

UiJ= l dl ! dll ^p^L (2) 

J J e (ri;r 2 )|ri - r 2 | 
where n^r-i) = \ipi(p, z, 0)| 2 is the particle density of QD, and eo(rx; r 2 ) denotes the position-dependent 
static dielectric constant. For the purpose of constructing approximate wave functions, we place the 
system in a large cylindrical confining box with the length L and radius R (L and R must be much 
larger than those of cylindrical Ge QDs). Here L = QOnm and R = AOnm are adopted. We solve the 
eigenfunctions of the effective-mass Hamiltonian by the Ritz variational method. The wave functions 
are expanded in a set of basis functions, which are chosen to be products of Bessel functions and sine 
waves ip n ,e,m(p, z,<f>) = Je(PnP)e li ^sin(k m (z + L/2)), where k m = mir/L,m = 1,2,3.., .In is the Bessel 
function of the order of i and f3 n R is the nth zero of J^. The same set of basis functions has been used 
by Marzin and Bastard[13] to calculate the quantum confined states in a conical QD. The expression 
of the matrix elements of the Hamiltonian of Eq. (1) can be readily obtained. Forty sine functions 
multiplied by fifteen Bessel function for each angular function (£ — or 1) are used to diagonalize 



the Hamiltonian. Fig. 2 show the lowest two energy levels of a single dot as a function of radius 
R for £ = 0. Note that the lowest energy level for £ = 1, not shown in Fig. 2, is higher than the 
energy level denoted by the dashed line. To determine physical parameters such as intradot Coulomb 
interactions, and interdot Coulomb interactions, we consider the dot A of radius R$ = 4.6nm and 
the dot B of Rq = AAnm. The ground state energies for the dot A and dot B are, respectively, 
O.lOSeV and 0.118eV. Fig. 3 shows the intradot Coulomb interactions and interdot interactions as 
a function of dot separation (ds), which is defined as the distance between the bottom of dot A and 
dot B. The solid line and dashed line denote the intradot Coulomb interactions of the dot B and dot 
A (Ub = ?>7.\-%meV and Ua = 36.437776^), the dotted line denotes the interdot Coulomb interactions 
(Uab)- Obviously, the intradot Coulomb interactions are almost independent of the dot separation. 
It is also found that the interdot electron hopping energy, t, is smaller than O.lmeV for ds > 16nm. 
This negligible t is due to high potential barrier between Ge QDs and Si0 2 - This implies that the 
probability of charge transfer from the dot A to the dot B is small. 



III. TUNNELLING CURRENT 

Now, we construct the following Hamiltonian to describe the system as shown in Fig. 1, 

H = Y e kdl,a,e a k,a,e + Y E i4,A,a + U i,j n i,v n 3,°' ( 3 ) 

k,cr,£ i=A,B;a i,j,a,a' 

k,<T,£,i 

where a) kai creates an electron of momentum k and spin a with energy e k in the £ electrode. d\ c 
creates an electron inside the SET's QDs with orbital energy E iy Uij describes the intradot and 
interdot Coulomb interactions and V4/,j describes the coupling between the band states and QDs. The 
Hamiltonian given by Eq. (3) is based on the Anderson model with two energy levels[14]. The 
time-independent tunnelling current from the left leads (source electrodes) can be expressed as (ref. 
[14]): 

— 9p r df 1 

J=—i:J ^r LiJ ( e )/m[-G< CT ( e ) + hAe)Gl a (e)}, (4) 
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where /^(e) = (exp[(e — H£,i)/k B T\ + is the Fermi distribution function of the I = L,R lead. 
r^j(e) = 2nJ2k I Vfc,i 1 2 ^(e — e k ) denotes the electron tunnelling rate from the ith dot to the I lead. 
Notations e and H denote the electron charge and Plank's constant. The lesser Green's function 
Gf a {e) and retarded Green's function G r ic {e) are the Fourier transformation of Gf a {t) = i(d\(0)di(t)) 
and G r ia (t) = -te(t)({d^{t),d\ a {0)}), where 9(t) is a step function, the curly brackets denote the 
anti-commutator, and the braket (...) represents the thermal average. Using Dyson equation, we 
obtain the lesser Green's function 

ImGf = _ r»(.)/,,.( e ) + r R1 (.)/ R1 (.) 2/m 

Substituting Eq. (5) into Eq. (4), we obtain the tunnelling current as 

J=Z 2e v fde r i|L (e)r i|fl (e) 

ft .^ B 7 2tt r i)L (e) + r i)fl (e) 1 ' 

For simplicity, these tunnelling rates will be assumed energy-and bias- independent. Therefore, 
the calculation of tunnelling current is entirely determined by the spectral function A = ImG r i a {e), 
which is the imaginary part of the retarded Green's function G r i a (e). The expression of retarded 
Green's function G r i a (e) can be obtained by the equation of motion technique[15-18]. The lowest 
order coupling between the electrodes and QDs will be considered in the calculation of Gg (e). The 
equation of motion for G r i a (t) leads to 

(e-E i + ^)G^( e ) = 1 + U^e) + U^G^e) + G r ij2 (e)), (7) 

where the two particle Green's functions G r u {e), G^i(e) and Gr[- 2 (e) arise from the particle correlation 
and satisfy 

(e - (E t + Ui) + *y)^(e) = N it . a + U tJ (G r m (e) + G\ ji2 {e)\ (8) 
(e - (Ei + Ui,) + i^G^e) = N ha + U t G r m (e) + U^G^e), (9) 

and 
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(e - (Ei + Uij) + i^)G r ij2 (e) = N jt . a + U t G r lJt2 (e) + U tJ G r w (e). (10) 



The notations Ni- a , Nj i<T and Nj_ a are the electron occupation numbers for the % dot. Note that 
% is not equal to j in Eqs. (7)-(10). Now the two-particle Green's functions are coupled to the 
three-particle Green's functions G r ijn (e) = (ni- a n j)r7 di i(7 ,dl a ), G r iji2 (e) = {ni-^Uj-adi^d]^), and 
G^- -(e) = (nj- a nj ta d ija , d\ a ). The equation of motion of the three particle Green's functions will lead 
to coupling with the four particle Green's functions, where the hierarchy terminates. Thus, these 
three particle Green's functions can be expressed in the following closed form 

1 " Ni-r. 



Gl jn (e) = N^N ha 



e-iEi + Ui + U^+i^ 

i Nj '~ a ) (a) 

e-(K + U i + 2U ii )+i^ h 



1-N<. 



G r iii2 (e) = Ni —tjNj _ CT ( — , -■ 



e-(E i + U j + 2U ij )+i i - 



+ Z ,„ , rr ^.rr i ■ tJ » ( 12 ) 



and 



1 - 



N- 

_i 1 ~ a ) (12) 

e-{E i + 2U ij + U i )+i^>' y J 

Substituting Eqs. (11), (12) and (13) into Eqs (9) and (10), after some algebras, the retarded 

Green's function G\ is given by 

(14) 

e-Ei- + e-Ei- 2U i:j + i^f 

+ N ._ r 1 -(% + jV ^)+% jV ^ 
e-Ei-Ui + i^ 

e- Ei-Ui-Uij e- Ei-Ui-2Uij 
Note that we ignored the electron-electron correlation effect on the tunnelling rates Tj in Eq. (14). 
The approximation used to obtain Eq. (14) is adequate for the Coulomb blockade but not for the 



Kondon effect. In Eq. (14) there are two main branches, one branch describes one-particle propagation 
process and the other describes two-particle propagation process. In each branch there are three 
channels with different probabilities. That the total probability of these channels is one indicates the 
satisfaction of sum rule 1 = j{de/n)\ImG r ic (e)\. Refs. [16] and [17] calculated the the tunnelling 
current of vertically coupled QD including interdot Coulomb repulsion and employed the mean-field 
scheme to solve the retarded Green's function, in which the resonant energy levels are found to be 
dependent on the occupation numbers of JVj and Nj. However, the physical interpretation for those 
resonant energy levels becomes ambiguity. In our calculation, the poles of Eq. (14) do not depend 
on electron occupation numbers, which determine the probability for each channel. The electron 
occupation numbers of Eq. (14) can be solved in a self-consistent way using 



N i>a is limited to the region of < JV i)(T < 1. Electron occupation number N A (N B ), according to Eq. 
(14), affects the tunnelling current through dot B (A). Eq. (14) is main result of this work. Based on 
Eq. (14), we calculate the current- voltage characteristic curves. 



In this section, we perform the detail numerical calculation of I-V characteristics. For simplicity, we 
assume that tunnelling rates of T 1 = 0.2meV and T 2 = 0.2meV are bias- independent. Meanwhile the 
chemical potentials of four electrodes are assumed to be 5meV below the energy level E A for zero bias. 
The electron Coulomb interactions can be regarded as bias- independent due to very strong confinement 
effect from the barrier height of V = 3.5eV. To determine the interdot Coulomb repulsion, we consider 
the interdot separation distance ds = 16nm, which offers Ua,b = 11.06meV. Substituting Eq. (14) 
into Eq. (15), the electron occupation numbers N A and N B (for fixed electron spin) could be solved 
in a self-consistent way. Figure 4 shows N A and N B as a function of applied voltage V a at zero 
temperature. We consider zero temperature through out this article. In order to exhibit the interdot 




(15) 



IV. RESULTS AND DISCUSSIONS 
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Coulomb repulsion effect, the dashed lines for Ua,b = are also plotted in Fig. 4. In the absence of 
Ua,b, the retarded Green's function of Eq. (14) has only two poles 



' i,—a 



1 — TV- N- 

Therefore, the staircase behavior of N is generated by the on-site Coulomb interactions £/;. Owing 
to T L — Tr, Ni reaches around 0.33 when the applied voltage crosses the energy level E i: but not 
sufficiently to overcome the charging energy C/j. This is the well-known Coulomb blockade effect. 
When the applied voltage overcomes the Coulomb blockade, iVj reaches to 0.5. These fractional 
occupation numbers are the typical statistic feature for an open system. 

For Ua,b 7^ 0, Ni shows a complicated behavior. We see that the "switch-on" voltage V a \ for 
B-SET is not changed in the presence of tunnelling current of A-SET. Nevertheless, the original 
occupation number N B = 0.33 (for Ua,b — 0) is reduced as the value of N B = 0.22. Once N B is finite, 
the probability of resonant energy level Ea is reduced. Consequently, the value of N A is suppressed. 
When the applied voltage V a reaches the resonant energy level E A + Uab, Na increases its value 
again. As the applied voltage is modulated to V a2 , the value of N B increases since the resonant level 
E B + Uab is opened. Na is slightly reduced at V a 2- The reduction of Na is small at V a 2 because the 
one particle process for Na is determined by the factor (1 — A^i_ CT )(l — N Bja N B _ a ) , where N Bja N B -„ 
is around 0.09. When the applied voltage is tuned to V a = A2mV for N A and V a = 53mV for N B , 
the two particle process is occurred for both Na and N B .The values of Na and N B that are larger 
than 0.33 are manifested results of exhibiting two particle process. Obviously, Na and N B are still 
suppressed by the interdot Coulomb repulsion , even though the applied voltages overcome the on-site 
Coulomb interactions. 

Once the electron occupation numbers Na and N B are determined, the tunnelling current can 
be calculated by substituting Eq. (14) into Eq. (6). Figure. 5 exhibits the tunnelling current and 
differential conductance (defined as dJ(V a )/dV a ) as a function of applied voltage V a . Diagrams (a) 
and (b) display the current through the dot A and dot B, respectively. We see that the interdot 
Coulomb repulsion not only creates new resonant channels for the source electrons transferring into 
the drain electrode, but also yields negative differential conductances. Results of Fig. 5 clearly exhibit 
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the suppression of tunnelling current by the interdot Coulomb repulsion. The negative differential 
conductance through the dot A(B) is generated by a fact that the applied voltage is not sufficient to 
overcome the Coulomb repulsion from the charge resided at the dot B (A). 

Finally, we study the dot B situated at the condition of T b ,l 3> T b ,r- From experimental point 
of view, this is due to extremely asymmetrical barrier width. For the dot B with the condition of 
T_B,L r_B,_R, the B-SET behaves as a closed system. We assume that the dot B can be manipulated 
as the following three charge states; the empty state \B,0), one electron state with spin up or down 
\B, t) (\B, I}) and two electron states \B, f , |). Figure. 6 shows the tunnelling current and differential 
conductance as a function of applied voltage for various charge states of the dot B; solid line denotes 
\B, 0), dashed line denotes \B, |) (\B, |)) and dotted line denotes \B, |, |). The separation of the first 
three peaks clearly displays the interdot Coulomb repulsion Ua,b- This indicates that the "switch-on" 
voltage of A-SET depends on the charge number of dot B. Obviously, the dot B situated in a closed 
system will significantly affects the "switch-on" voltage of A-SET. This also demonstrates that an 
SET is a highly sensitive charge detector. 

V. SUMMARY 

For parallel double SETs implemented by QDs embedded into Si0 2 , the interdot Coulomb repul- 
sion is more important than the interdot electron hopping as a result of high potential barrier. The 
interdot Coulomb repulsion not only creates new resonant energy levels but also produces negative 
differential conductances. Nevertheless, note that the "switch-on" voltage of SETs is not changed 
when both two SETs behave as an open system. If one of SETs such as B-SET in this study behaves 
as a closed system, and then the role of B-SET acts like a trap. The "switch-on" voltage of A-SET 
will be dependent on the charges resided at the QD of B-SET. This implies that the precisely control 
of barrier width in the implementation of circuit integrated by SETs is very crucial. 

This work was supported by National Science Council of Republic of China Contract No. NSC- 
93-2215-E-008-014 and NSC 93-2120-M-008-002 
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Figure Captions 

Fig. 1: Schematic diagram for parallel double single electron transistors. 

Fig. 2: The lowest two energy levels of a cylindrical Ge/Si0 2 QD for £ = as a function of dot 
radius i?o- 

Fig. 3: Particle Coulomb interaction strengths as a function of dot separation ds. 
Fig. 4: Electron occupation numbers as a function of applied voltage. 
Fig. 5: Tunnelling current and differential conductance as a function of applied voltage. 
Fig. 6: Tunnelling current and differential conductance through the dot A for various charge states 
of the dot B. 
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